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Abstract 



The composite operator effective potential is compared with the conventional 
Dyson-Schwingcr method as a calculational tool for (2+l)-dimensional quantum 
electrodynamics. It is found that when the fermion propagator ansatz is put directly 
into the effective potential, it reproduces exactly the usual gap equations derived 
in the Dyson-Schwinger approach. 



1 Introduction. 



Studies of iV-flavour Quantum Electrodynamics in three dimensions (QED 3 ) have 
long been made in order to learn about dynamical mass generation (and hence 
chiral symmetry breaking) in a relatively simple context jl], [^, ||, Initially the 
hope was to gain insight into dynamical symmetry breaking in other asymptotically 
free theories, in particular Yang-Mills theories in four dimensions. Subsequently 
there has been interest in QED3 as a model for high-T c superconductors ||. Here 
the broken, massive phase of QED3 gives a mechanism for spontaneously breaking 
the £7(1) gauge symmetry which is necessary for superconductivity, and deviations 
from the trivial infra red fixed point could be responsible for the observed non- Fermi 
liquid behaviour in the normal phase of these materials || . 

Considering N fermion flavours allows one to perform a large N expansion and 
keep only dominant terms in this limit. This is a convenient calculational tool for 
it admits a form for the full gauge boson propagator (namely that with massless 
fermion loops summed Q) without requiring that the Dyson-Schwinger equation 
for the propagator be solved. There is a further significance of working with iV 
flavours in the possibility that there exists a critical number of flavours f?| above 
which no symmetry breaking occurs. Much of the evidence for the existence of a 
critical flavour number lies in numerical work and the analytical studies that have 
been made require great simplifications to be tractable; indeed there is evidence for 
no critical behaviour in the flavour number (see e.g. Q). 

Non-perturbative treatments of QED3 are normally based on solving systems 
of Dyson-Schwinger equations. To make the systems soluble the equations for 
the gauge boson self energy are dropped in favour of the propagator resummed in 
large N; the equations for the vertex are replaced by an ansatz discussed further 
below. This leaves only the Dyson-Schwinger equation for the fermion propagator 
to be solved; the approach taken is to substitute a (rather general) ansatz for the 
propagator into the equations and then solve for the functions appearing in the 
ansatz. 

The system of Dyson-Schwinger equations for a theory can in principle be de- 
rived from a composite operator effective action |^| via a well-defined minimization 
technique. It is the purpose of this paper to consider whether the same physics is 
found if one puts the ansatz directly into the effective action. The same minimiza- 
tion technique can then be used to derive the "gap equations" for the functions 
in the ansatz directly from the effective action, and these can then be compared 
with the conventional Dyson-Schwinger results. The system of Dyson-Schwinger 
equations used is heavily truncated and the effective action is computed to two loop 
level only; it is therefore not immediately obvious that the results of the two compu- 
tations will coincide. Indeed in general the calculations may give different answers 
(e.g. j^]), and comparing the methods can allow one to estimate the limitations 
of the approximations used. Since the approximations described above are used 
almost universally in work on QED3 it is important to compare the two methods 
within the those limits. It will be shown that the composite operator effective action 
method gives exactly the usual (Dyson-Schwinger) gap equations, after simplifica- 
tion of the "raw" equations which come directly from the minimization procedure. 
The amount of computational effort required is rather surprising, and may have 
important consequences for future work at finite temperature; it is probable that 
finite temperature Dyson-Schwinger approaches would be more economical. 

This paper is organized as follows: in section^, the basic formalism for composite 
operator effective actions is reviewed, rewritten for fermionic fields; section |3| con- 
tains a comparison of the effective action and Dyson-Schwinger equation approaches 
and in section |J are the details of the effective action calculation. Conclusions are 
appended to these sections. 
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2 The Composite Operator Effective Action. 

The composite operator effective action || is a generalization of the conventional 
(quantum) effective action; for fermionic fields it depends not only on ip{x) and ip(x) 
— possible expectation values of the quantum fields ^(x) and ^f(x) — but also on 
S(x,y), a possible expectation value of the time ordered product T {'5 (x) 1 ^ (y)} . 
Physical solutions require that the effective action be stationary with respect to 
variations in the expectation values: 



Sip(x) 
5ip(x) 



0, 
0, 



(1) 



5S(x,y) 

For this reason the formalism is especially suitable for studying dynamical symmetry 
breaking, for although the first two of equations ([!]) may only have the symmetric 
solution ijj(x) = ijj(x) — 0, symmetry breaking solutions may exist for the last 
equation. 

First we define the generating functional Z[r],rj, A] (J = J d 3 x): 
Z[v,rj,A] = exp{iW[r),rj,A]} 

+ / (rj(x)^f(x) + l>(x)ri(x)) 



D^D^ exp < i 

I L J x 

nx)A(x,y)^(y)]\ , (2) 

where {77 , 77} and A are fermionic and bosonic source currents respectively. The '5 
and W integrations are functional, and I[^f, Vt] is the classical effective action. The 
composite operator effective action is the double Legendre transform of W[rj,fj, A]. 
Defining 



SW[ri,T], A] 

5rj(x) 
5W[r],rj, A] 

8rj(x) 
6W[r],rj, A] 



= ^(x)^(y) + S(x,y), (3) 



SA(x,y) 

one can eliminate 77, rj and A in favour of tp, tp and S in the Legendre transform 
T[^,S] = W[t},%A]- [ @(x)ri(x)+Jj(x)il>(x)) 



i/j(x)A(x,y)xp(y) - S(x,y)A(y,x). (4) 

i x J y J x J y 



It is then easy to show that 

sr[y>,^,s] 

6ip(x) 

5%j)(x) 

ST[^,S] 
6S(x,y) 



rj(x) + / ip(y)A(y,x), 

-r){x) - / A(x,y)ij}(y), 
Jy 

-A(x,y). (5) 
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Physical processes correspond to vanishing sources, so equations (^) provide a 
derivation of the stationary requirements, equations (|l|). 

Since we will be interested only in translation-invariant solutions we can set ip 
and tp to constants and take S(x,y) to be a function only of the separation x — y. 
This leads to a generalization of the effective potential || , defined by 

V[^,S\ f = -T[^,S]. (6) 

J X 

Then a series expansion for V[ip, tp,S] is relatively easy to find ||; the specific 
expansion for QED3 is given in section |3|. 

Composite operator effective actions depending on possible expectation values 
of irreducible three-point, four-point etc. functions can be defined simply by taking 
further Legendre transforms, however the "bilocal" effective action defined above is 
sufficient to study QED3. 

3 The Effective Potential and Dyson— Schwinger 
Equations in QED3. 



The Lagrangian density for iV-flavour QED3 is as follows: 

JV 

F 2 -4- i ' 
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1 N - 

C = —Fl v +iY,^kl»V^ k + C GFl (7) 



fe=i 



where and are the usual U(l) curvature tensor and covariant derivative 
respectively; the fermions are represented by four-component Dirac spinors so 
that the theory has a chiral symmetry and can exhibit mass generation gj. The 
gauge fixing term is not written explicitly: throughout the calculation we shall 
work in Landau gauge. The series expansion for the composite operator effective 
potential in momentum space is then || (J, = J d 3 k/(2Tr) 3 ): 

V[S,D] = iN fix {HS- 1 (p)S (p)} + S^(p)S(p) - 1} 
Jp 

- % - J tr {]n[D- 1 (p)D (p)] + D 1 {p)D{p) - 1} 

+ V 2 [S,D]. (8) 



The functions S and D are candidate full non-perturbative propagators for the 
fermions and gauge fields respectively, and the subscripts denote their bare coun- 
terparts. The effective potential is minimized with respect to (functional) variations 
in S and D to determine the physical non-perturbative propagators for the theory, 
as described in section 0. — V^S", D] is the sum of all two-particle irreducible vac- 
uum graphs with the propagators set equal to S and D and with bare (not dressed) 
vertices. To consider more general vertices one would have to got to a "trilocal" 
effective action (i.e. three Legendre transforms), and require that this new object 
be stationary also with respect to variations in a vertex ansatz. 

The graphs contributing to V2 [S, D] up to three loop level are depicted in figure 

There are also two-particle reducible graphs which would have to be considered 
in an ordinary (one-particle irreducible) effective action calculation but which do 
not contribute to V% [S, D] ; these are depicted up to three loop level in figure |^. 

We will truncate the series at the two loop level, as this will be sufficient for 
comparison with the usual Dyson-Schwinger method for QED3 . The two-loop term 
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Figure 1: Two-particle irreducible graphs contributing to V2, up to three loop level. 




Figure 2: Two-particle reducible graphs not contributing to V2 (to three loop level), 
but which would contribute to the one-particle irreducible effective action. 



in V2 follows: 

e 2 N f f 

V 2 [S,D] = —J J^T{^S(p)T u S( P + k)D^(k)}. (9) 

For the purposes of comparison with the Dyson-Schwinger method we keep the 
"general" vertex T; when we turn to the effective potential calculation proper in 
section ^ this will be reset to the bare vertex. On demanding that V[S, D] be sta- 
tionary with respect to variations in S and D the usual Dyson-Schwinger equations 
are obtained: 

S-'ip) = So X (p) +e 2 J^ {S(k)^D^(q - k)r„} ; (10) 
D-l{p) = D^\ v {p)-e 2 N f {^S{k)Y v S{k+p)}. (11) 

J k 

Now these equations can be solved for S and D (with some ansatz or a further 
Dyson-Schwinger equation for the vertex) to give the full non-perturbative prop- 
agators which minimize the effective potential. Since the interest is in possible 
dynamical mass generation in the fermion sector we use an ansatz for the fermion 
propagator and solve for the functions appearing therein: 

S- 1 (p) = -i(A(p) j+B(p)). (12) 

Here A and B are both assumed to be scalar functions, so that this ansatz is rather 
general in that it only omits possible parity-breaking mass terms. 

The gauge propagator is parameterized by the vacuum polarization, and we 
approximate this function by the leading order truncation in 1/N Q: 

D *M—?T*r r { , -- t ?)-' (13) 
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2 N 



The only requirement put on the vertex function is that it is dependent on the 
magnitudes alone of the momenta at the vertex, but it must also be consistent with 
the Ward-Takahashi identity pCj : 

T fi {p,k) = -i lli G(p 2 ,k 2 ). (14) 

The full vertex ansatz (|lj) appears only once in V2 to avoid double-counting. This 
treatment of the vertex is consistent with that in Dyson-Schwinger calculations; 
when we turn to the effective potential calculation, the function G will be set to 
unity. In order to take account of higher loop effects in the vertex, one has to 
consider the "trilocal" effective action discussed briefly above. 

When these functions are substituted into the Dyson-Schwinger equations ( |l0| , 
pr] ) the following gap equations for A and B can be derived: 

" 1 f°° „ kACk^Giq 2 ^ 2 )^ „ 



B(q) 



a 1 



tt N q 



dk 



kB(k)G(q 2 ,k 2 
KJJi) 



with the integrals given by 
Ii(k 



,q;a) = 



a 



In 



-\k 2 -q 2 \(k 
a 

1 / , o o\ 



a (k -\ 

q-\k-q\) 
k 



In 



q - 



l2(k,q;a) = 4 In 



k + q 



and 



\k-q\ 



B 2 (k). 



-T 2 (k,q;a); 

■q-\k-q\ 
+ 2kq 
-In 



(15) 



k — q 



(16) 



(17) 



JC(k) = A 2 (k)k 2 

These equations can now be solved (using analytical and numerical methods) for A 
and B to yield the full fermion propagator which minimizes the effective potential. 

The approach taken in this paper will differ from the conventional Dyson- 
Schwinger method outlined above. We shall use the ansatz ( |l^ ) directly in the 
effective potential, equation (||),and demand that V[A, B] be stationary with re- 
spect to (functional) variations in A and B. We can then derive the "gap equations" 
for A and B and compare them with the equations found in the Dyson-Schwinger 
approach, equations (|l5|). 



4 Gap Equations from the Effective Potential. 

The effective potential comprises three natural parts; the trace over the fermion 
propagator terms, the trace over the gauge boson propagator terms, and the sum of 
two-particle irreducible vacuum graphs. The second of these will play no part in the 
variation to be described: we will again use the leading order (in 1/-/V) truncation 
for the gauge propagator (|l^). After substituting (|l2|), the variation of the pure 
fermion part is given by: 

SV f [A,B] = ANi A 2 {q)q i {A{q) - 1) + B 2 {q)q 2 {A{q) + 1) 

5A(q) ~ (2tt)3 K?{q) 

SVf[A,B] _ 4Ni A(q)B( q ) q 2 (A(q)~2) + B 3 (q) 

SB(q) (2tt) 3 JC 2 (q) ' 1 ' 
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When these variations alone are set to vanish, the equations have the solutions 
A = 1, B = 0, consistent with the free propagators. The variation of the loop part 
is: 



SV 2 [A,B] 
SA(q) 



SV 2 [A,B] 
SB(q) 



ie 2 N2A{q)B{q)q 2 f B(k) 



1 



2tt 3 
ie 2 N 
~2^~ 
ie 2 N 

2tt 3 

ie 2 N 2A{q)B{q) f A(k) 



JC 2 (q) 
2A 2 (q)q 2 
IC 2 (q) 

2B 2 (q) 

JC 2 (q) JC(q) 



k K{k) {k-q) 2 l + a/\k-q\ 



JC{q) 
1 



A(k) 



q ■ (k - q)k ■ (k - q) 



k K{k) (k - g)4 
B(k) 2 



1 + a/\k - 



k K{k) (k-q) 2 l + a/\k-q\ 



q ■ (fc - q)k ■ (k - q) 



2tt 3 K 2 (q) J k JC{k) (k - q) 4 1 + a/\k - q\ 



(19) 



Computing the resulting angular integrations, performing some simplifications, 
and setting the variations to vanish we obtain the "raw" gap equations: 







A 2 (q)q i (A(q)-l)+B 2 (q)q 2 (A(q) + l) 



2a 

Nn 2 ' 



2a 



Nir 2 

= A(q)B(q)q 2 (A(q) 
2 a 



(A 2 (q)q 2 -B 2 (q)) \ dk J 2 (k, q; a), 

2a 



2A(q)B(q)q 2 \ dk J x {k, q; a) 
(20) 



B 3 (q) + ^(A 2 ( q )q 2 -B 2 (q)) J dkj^q-a) 



2A(q)B(q) J dkj 2 (k,q;a) 



Nir 2 ' 

where the integrals are given by: 



(21) 



Ji(k,q,a) 
J 2 (k,q;a) 



k2 B(k) 2 



hi 



2k 



K{k) kq 
2 A(k) f(k 



k + q + a 
\k — q\ + a 

2_ q2)2 _ a 4 



In 



K.(k) V 

(k 2 -q 2 ^ 2 

Akqa 2 
(k 2 - g 2 ) 2 r 
Akqa 



Akqa 2 
k + q~ 



In 



\k-q\ 



\k-q\ 
1 



W q [k + q -\ k - q \ ] -2 



\k — q\ k + q 



(22) 



Since the momentum integrals appearing in equations (|2C|, 21) are the same we 
can eliminate each of them in turn: after a little simplification this yields the usual 
gap equations obtained from the Dyson-Schwinger approach, equations ( p"5| ) , up to 
the replacement G — > 1 which is imposed by considering only a "bilocal" effective 
action. 



5 Concluding Remarks. 

The fact that the composite operator effective potential yields the same gap equa- 
tions ( |l5|) whether the ansatze are put in before or after the functional variation 
demonstrates explicitly that it is equivalent to the Dyson-Schwinger method for 
QED3 . It is by no means obvious from the "raw" gap equations ( f20| , pTl) that they 
will reduce on simplification to the normal Dyson-Schwinger results; since the effec- 
tive potential and Dyson-Schwinger methods are frequently applied to QED3 it is 
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important to have checked their equivalence within the framework of the approxima- 
tions normally used. The consistency here is a basis for confidence that the ansiitze 
used are sufficiently consistent and general: in other work || it has been found 
that overly restrictive ansiitze lead to differences between Dyson-Schwinger type 
approaches and the results of composite operator effective action computations. 

To better model high-T c superconductivity, finite temperature effects should 
really be included. The conclusion that can be drawn from the present calculation 
is that using the composite operator effective potential for such work (see e.g. pJ|) is 
not likely to yield additional information nor to simplify the computation compared 
with a finite temperature Dyson-Schwinger method. 

At both finite and zero temperature the choice of vertex function remains the 
major problem with this sort of calculation: how to make an ansatz that is consistent 
with the Ward-Takahashi identity, when that identity requires the full fermion 
propagator. This problem is not resolved by using the "bilocal" composite operator 
effective potential; however, in principle this problem could be solved completely 
by using a "trilocal" composite operator effective action, leading to a prediction 
for the vertex. However the calculation would be extremely difficult. The most 
promising way to attack this particular problem is by using the non-local gauge , 
in which the wave-function renormalization A is forced to be unity. This approach 
does not work well within the composite operator effective potential formalism, for 
the resulting analogue of equation (|2f]) for B is intractable. 
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